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Abstract 



We compute bulk properties of Heisenberg spin- 1/2 ladders using Rayleigh- 
Schrodinger perturbation theory in the rung and plaquette bases. We formu- 
late a method to extract high-order perturbative coefficients in the bulk limit 
from solutions for relatively small finite clusters. For example, a perturbative 
calculation for an isotropic 2 x 12 ladder yields an eleventh-order estimate of 
the ground-state energy per site that is within 0.02% of the density-matrix- 
renormalization-group (DMRG) value. Moreover, the method also enables a 
reliable estimate of the radius of convergence of the perturbative expansion. 
We find that for the rung basis the radius of convergence is A c ~ 0.8, with A 
defining the ratio between the coupling along the chain relative to the coupling 
across the chain. In contrast, for the plaquette basis we estimate a radius of 
convergence of A c ~ 1.25. Thus, we conclude that the plaquette basis offers 
the only currently available perturbative approach which can provide a re- 
liable treatment of the physically interesting case of isotropic (A = 1) spin 
ladders. We illustrate our methods by computing perturbative coefficients 
for the ground-state energy per site, the gap, and the one-magnon dispersion 
relation. 



I. INTRODUCTION 



Heisenberg spin ladders have been extensively studied via a number of powerful theo- 
retical methods including direct diagonalization JI], quantum Monte Carlo [IJ and density 
matrix renormalization group (DMRG) @. Such investigations have been complemented by 
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development of simpler, often analytic, treatments whose role is to help elucidate the physics 
underlying the numerical results. These approaches are typically based on perturbation the- 
ory and almost invariably assume a "dimer" or "rung" basis (see Fig. [1]). Note that in the 
rung basis, the two-leg-ladder Hamiltonian becomes diagonal in the limit in which the cou- 
pling along the chain (Jy) becomes negligible relative to the coupling across the chain (J±). 
In this basis, the intra-chain couplings are, thus, treated perturbatively. Such perturbative 
expansions are routinely assumed to be meaningful, even in the region where the inter-chain 
and intra-chain couplings are comparable (i.e., J\\/J± ~ 1) 

In the present work we present a perturbative treatment of two-leg spin ladders imple- 
mented using a "plaquette" basis which, as we have argued elsewhere [f|f|, is particularly 
well-suited to the study of such systems for a variety of reasons. Our method exploits the 
simplicity of these quasi one-dimensional systems to straightforwardly extract perturbative 
information about the bulk system from ordinary Rayleigh-Schrodinger perturbation theory. 
In so doing we avoid the more general, and, thus necessarily more complicated, diagrammatic 
methods of "connected graph expansions" 0. Specifically we can, with modest computa- 
tional effort, determine the bulk values of the ground state energy per site, the spin gap, 
and the one magnon dispersion relation to eleventh, fourth and third order, respectively. 
Moreover, by performing similar calculations in the rung basis, we can make reliable quanti- 
tative comparisons of the radii of convergence of the perturbative expansion in the two bases. 
Indeed, we determine that the radius of convergence in the rung basis is A c = J\\/J± — 0.8 
while for the plaquette basis, where, as will be discussed below, A scales the inter-plaquette 
coupling, we find A c ~ 1.25 to 1.4 . Clearly, the plaquette basis is strongly preferred for the 
treatment of the physically interesting || case of isotropic (i.e., A = 1) spin ladders. 

Our paper is organized as follows. In Sec. |I| we briefly describe the rung and plaquette 
bases and outline how to apply Rayleigh-Schrodinger perturbation theory to finite spin 
ladders in these bases. In Sec. [H| our method for extracting perturbative information for 
bulk systems from the finite-ladder results is outlined. Our numerical results appear in 
Sec. |V| where we also address the issue of the radii of convergence in the rung and plaquette 
bases. The paper concludes with a summary section. 

II. PERTURBATION THEORY IN COUPLED BASES: FINITE SYSTEMS 

A. The Rung Basis 

For a two-leg ladder, the dynamics are described by the Heisenberg Hamiltonian with 
nearest-neighbor antiferromagnetic coupling: 

H=J\\Y,Si-S j + J±Y, S i- S 3> C 1 ) 

~ I 

where J±(J\\) is the strength of the antiferromagnetic coupling along the rungs(chains) of 
the ladder |TJ]. In the strong-coupling limit (J\\/J± — 0) rungs interact very weakly with 
each other and the above separation of the Hamiltonian is a natural one. Indeed, by fixing 
the large scale in the problem to be one (J± = 1) the Hamiltonian can be separated into 
the sum of Hq, which contains the intra- rung interactions, and a potential V — scaled by the 
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small parameter A = J\\/J±, responsible for inducing interactions between rungs. In this 
way, the dominant intra-rung Hamiltonian can be written as (see Fig. ||) 

N R 

H = J2hP; = S?> • S<"\ (2) 

r=l 

while the inter-rung potential (XV) is expressible in terms of 



v = J2 [s! r) ■ sS r+1) + s< r) ■ sf +1) ] . (3) 



r=l 

Note that r labels a specific rung and the sum runs over all Nr rungs in the ladder. 

The rung basis is particularly suitable in the strong-coupling regime, as H is diagonal in 
this basis. Indeed, eigenstates(eigenenergies) of H are a direct product (sum) of single-rung 
states (energies) of the form: 

\<f>tm) = I (sis 2 ) lm) ; q = li(£ + 1) - 3/4 . (4) 

Here the two spins along the rung are coupled to a total angular momentum £, which can 
be zero or one, and projection m. What remains to be done is to evaluate matrix elements 
of the potential in the rung basis. To do so, we rely on the same angular-momentum 
techniques fl,|10|l that were used to evaluate matrix elements of the Hamiltonian in the 
plaquette basis plH. We obtain 



(5) 



(f (1 ¥ (2) jm|y(l,2)|£ (1 ¥ 2) jm> = ±3£( 1 )£( 2 )£'( 1 )£'( 2 ) 
£W f(D 1W£(D f(D 1 W£(2) f(2) i 
f ( 2) £ ( 2 ) j}\i/2 1/2 1/2 J 1 1/2 1/2 1/2 



Note that we have defined x = \J1x + 1 and that the +/— sign in the above expression should 
be adopted whenever £^ + + j is even/odd. Moreover, the above expression vanishes 
whenever + £^ + + is odd. Aside from simple phases and numerical factors, 
the matrix elements of V(l, 2) depend on scalar functions (i.e., independent of m) known as 
Racah coefficients; here we cast the matrix elements in terms of the more symmetric 6— j 
symbols PJlO . 



B. The Plaquette Basis 

For an isotropic (J\\ = J± = 1) two-leg ladder, an efficient decomposition can be carried 
out in terms of distinct pairs of adjacent rungs; these are the plaquettes. In such a decom- 
position, a 2 x 4 ladder, for example, would be viewed as a pair of interacting plaquettes. 
The ladder Hamiltonian can be written as the sum of Hq, which contains the intra-plaquette 
interactions, and V, which includes the interactions between plaquettes. The intra-plaquette 
Hamiltonian H is expressible as (see Fig. [I]) 



Hn 



N P 

E 

P =i 



h 



(p) 



h 



si' 



(p) 



+ sl p) 



(6) 
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where p labels a specific plaquette and the sum runs over all Np plaquettes in the lad- 
der. As was shown in Ref. ||, H is diagonal in the plaquette basis. The single-plaquette 
wavefunctions and energies are of the form: 

\(/> a ) = |(s 1 s 4 )£ 14 , O2S3K23; 3 m ) , (7) 
e a = + 1) - iu{ii4 + 1) - + 1)] , (8) 

where a= {iu, £23, j, m} denotes a generic quantum label. In the plaquette basis the two 
diagonal pairs of spins are coupled to well-defined total angular momenta, £14 and £ 23 , which 
can equal zero or one. These two link angular momenta are in turn coupled to a total 
plaquette angular momentum j with projection m. In this scheme, the interaction V is 
given by 

Np 

v = J2 [s 3 p) ■ s? +1) + sf • s ( 2 p+1] ] ■ (9) 
p=l 

Matrix elements of V in the plaquette basis are computed in terms of five Racah coefficients, 
as was shown in Ref [01. 



C. Rayleigh-Schrodinger Perturbation Theory 

In this section we review some well-known results in perturbation theory |Tl,[l2l. To for- 
mulate Rayleigh-Schrodinger perturbation theory for finite ladders, we let the Hamiltonian 
H — > H(X) = H + AV^ and make a consistent expansion of the energy and wavefunction in 
powers of the "small" parameter A. That is, the exact wavefunction of H(X) is expanded in 
a power series in terms of the eigenstates of Hq 



l*«> = £A B £48l**> 



(10) 



n=0 



while the energy is expanded as 



(11) 



n=0 



Note that A^ = 5 n o, A^ = S aj 3, and Ef® = e a is an eigenvalue of the unperturbed Hamilto- 
nian Ho. It is straightforward to derive a recursion relation for the above coefficients which 
is well-suited for implementation on a computer to arbitrary order in perturbation theory. 
We obtain 



An) 



n-l 



m=l 



(12a) 
(12b) 



These formulas are very simple to apply to the ground state of a 2 x Nr ladder consisting 
of Np = N R /2 plaquettes. The unperturbed state is given by 
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$J= ) = |0O)®|0O)-..|0O) 



(13) 



In the rung basis |0o) represents the singlet state, while in the plaquette basis it is the 
1^14 = ^23 = 1, j = 0) single-plaquette wavefunction. The situation becomes slightly more 
complicated for the (J = 1) one-magnon states. At zeroth order this state is A-fold degen- 
erate (with N = Nr or N = Np) so Eq. [E| is no longer applicable. However, if we assume 
periodic boundary conditions, the total linear momentum becomes a good quantum number 
and serves to distinguish among the N degenerate states, which we now write as: 

»&) . (14) 

where the linear-momentum label runs from u — 0, . . . , N — 1 and 

|$S=i>H0o)®|0o>..-|0i>...|0o>- (15) 

Now |0i) is the triplet state for the jth rung, or the \£i4 = h3^= j = 1) single-plaquette wave 
function for the jth plaquette. In this way, Eq. [12| can be used for each individual state in 
the degenerate band, as states with different linear momenta cannot mix at any order in 
perturbation theory. We have systematically determined that, for sufficiently small values of 
A, our perturbative calculations agree essentially perfectly at large order with exact results 
obtained by direct diagonalization or Lanczos techniques. 



Hexp 



III. PERTURBATION THEORY IN COUPLED BASES: BULK SYSTEMS 

The essential concept that we use to extract perturbative quantities in the bulk limit 
from calculations for finite systems is that of "connectedness". It is well known in many- 



body theory (see, e.g., Refs. p| JTI|jr3]JI^| ) that any given perturbative contribution to the 
coefficients (-Ej= ) appearing in the expansion for the ground-state energy [see Eqs. (0,0)] 
consists only of spatially connected terms. This is also the basis of the connected-graph 
expansions alluded to in the Introduction. As is well-known flIT] , |i3| , |T4|| , both connected 
and disconnected terms appear in Rayleigh-Schrodinger perturbation theory. However, the 
disconnected terms are canceled by contributions included in the subtractions appearing 
in the evaluation of the A^I, as shown in Eq. ([12]). The linked-cluster expansions were, 
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of course, developed to avoid dealing with disconnected terms which would eventually be 
canceled anyway, the price being much greater conceptual complexity. 

In this work, we employ Rayleigh-Schrodinger perturbation theory and find that the 
computational penalty for not excluding disconnected contributions from the outset is in- 
significant compared to the advantages arising form the simplicity of the formulation. Still, 
we must exploit the connectedness of the surviving contributions to extract bulk properties. 
In the present context, connected terms are those for which a single, finite cluster of adja- 
cent plaquettes can be identified such that: 1) all interactions involve plaquettes within that 
cluster and 2) all plaquettes in the cluster are involved in one or more interactions (note 
that in the present discussion we make explicit reference to plaquettes; however, all the 
arguments remain valid in the rung basis, unless stated otherwise). It is apparent that, at 
order n in perturbation theory, the largest such cluster contains n + 1 plaquettes. Hence, we 



5 



may conclude that, at nth order, no new kinds of clusters appear as we change the number 
of plaquettes, N, provided N > n + 1. Instead, only the number of times that each type of 
cluster appears changes with N. As the energy must scale at most with N, we may conclude 
that the ground-state energy is generally expressible as: 

E { ]2 (N)=A + B N, (16) 

where Aq and Bq are "intensive" coefficients, which are independent of N, for N > n + 1. 
This simple, physically reasonable behavior — and more specifically the absence of terms 
depending on N 2 , N 3 , . . ., is again related to connectedness |Tl|JT^] . Indeed, it was in order 
to eliminate terms with "anomalous" iV-dependence that the linked-cluster expansion was 
first proposed. 

It is clear that, for the case of periodic boundary conditions, we may set A = 0. Thus, to 
nth order in perturbation theory, the ground-state energy-per-site of an iV-plaquette ladder 
becomes independent of the number of plaquettes for N > n + 1. Obviously then, the nth 
order energy-per-site in the bulk limit can be expressed as: 

^go/site = 1 -E^ (N = n + l)/(n+ 1) , (17) 

where Ej2 (N) is the singlet ground state for an iV-plaquette ladder, assuming periodic 
boundary conditions. 

The structure of the unperturbed one-magnon states (l^j^)) or, equivalently, the fact 
that the gap is an intensive quantity, ensures that, in the bulk limit, the energy-per-site 
equals that of the ground state. Hence, we may conclude that all one-magnon states have 
the same extensive contribution as the ground state. Therefore, 

E i J l 2 1 (v,N) = A^ + B N. (18) 

Again, to nth order, A^ can be determined from studying a finite ladder with at least 
N = n + l plaquettes. For u — 0, v4 ( ^ =0 ' ) is just the spin gap at order n and its value in the 
bulk limit is found simply by taking the difference between the nth order energies for the 
lowest one-magnon state and the ground state where both are evaluated for an N = (n + 1)- 
plaquette system with periodic boundary conditions: 

A (n) = E i ]2 1 (is = 0,N = n+l) - E ( j2 (N = n + l) . (19) 

We use this formula to determine the perturbative contributions to the spin gap in the 
bulk limit presented in the following sections. Note that in the rung basis, it is the linear 
momentum corresponding to u = N/2 (i.e., k = n) which generates the lowest one-magnon 
state. In this case one must employ a finite ladder containing an even number of rungs. 

It is also easy to show that the nth order one-magnon dispersion relation has the following 
form: 

n 

A (n) (£;) = ]T A( n) cos(fcr) , (20) 

x=0 

where < k < 27r is the linear momentum in lattice units. Note that that in the plaquette 
basis the gap is given by A^ n ^ = A^ n ^(fc = 0) while, in the rung basis, A^ = A^ n \k — it). It 
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might be anticipated that, at order n, the n + 1 coefficients which specify the full dispersion 
relation at this order could be extracted directly from the values of the differences 

A (n) (A; = 2wk/N) = Efl x {v, N = n + l) - E^N = n + l) . (21) 

However, a moment's thought reveals that some of these one-magnon states are degenerate 
and, thus, energies for systems larger than N = n + 1 are typically required. Indeed, for a 
2x8 (eight-rung) system, a perturbative calculation in the rung basis determines the bulk 
values of the energy-per-site and the gap up to 7th-order in perturbation theory, but only 
to 4th-order for the dispersion relation. Nevertheless, it is always possible to ascertain the 
width of the one-magnon band at n-th order from calculations for systems with N = n + 1 
plaquettes, provided N is even. 

Finally, it should be noted that the structure of the ground state in the plaquette basis 
(and, indeed, in any basis in which the unperturbed ground state is made out of singlets) 
affords a less obvious but more efficient means of extracting the coefficients that determine 
the ground-state energy in the bulk limit. To see this, we first note that |$j=o) contains only 
j = plaquettes [see Eq. (|I3"])1. Since the inter-plaquette interaction V, specified in Eq. (|J), 
consists of rank-1 operators in the space of a single plaquette, any plaquette appearing in a 
connected cluster must experience at least two interactions in order that it be returned to 
its original j = state. For this reason, it can be shown that EjZq(N), when computed with 
open boundary conditions, may be expressed as: 

[n/2] 

E { ;2 (N) = J2(N-m)W^(N = m + l) , (22) 

m=l 

where [n/2] is the largest integer less than or equal ton/2. Then, for example, the quantities 
Wi 5) and Wi 5) can be extracted from fifth-order calculations for systems with only two and 
three plaquettes. (Note that, in contrast, the extraction of fifth-order coefficients (Ej= ) for 
the ground-state energy using periodic-boundary conditions and Eq. flTTj ) requires calcula- 
tions for systems with six plaquettes.) Then, the nth-order energy-per-site in the bulk limit 
is easily found to be 

1 [n/2] 

Ego/site =-£ W^(N=m+l) . (23) 

4 m=l 

We use this formula to determine the bulk perturbative ground state energy contributions 
presented in the following sections. 



IV. RESULTS 
A. The Rung Basis 

We first demonstrate our method for extracting perturbative energies in the bulk limit 
from Rayleigh-Schrodinger calculations for finite systems. Fig. ^| shows rung-basis calcula- 
tions for the ground state energy and for the spin gap for systems with 4, 6, and 8 rungs. 



These calculations were done for A = 0.25 so convergence to DMRG values [15] is very fast 
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Fig. ^ shows that, at orders n = 2 and 3, all calculations give the same per-rung energy 
and spin-gap contributions, as we would expect from the arguments of the preceding sec- 
tion. At n = 4, the 6- and 8-rung values coincide, but now differ from the 4-rung results, 
again consistent with our expectations. Finally, the 6- and 8-rung values begin to differ 
for n > 6. Note that our 7th-order perturbative result for the spin gap agrees with the 
DMRG value to about 0.07% (for the energy per site both results are indistinguishable to 
five significant figures). The left-hand side of Fig. |2] also demonstrates that, when Eq. ( |2"3"D is 
used, the contributions to the ground-state energy extracted from 4-rung calculations with 
open boundary conditions agree with those obtained from an 8-rung ladder with periodic 
boundary conditions through n = 7, as expected. This result embodies the power of the 
method, and the high degree of symmetry of the unperturbed ground state and the resid- 
ual interaction. Indeed, an 8-rung calculation affords the computation of the ground-state 
energy — up to 15th order in perturbation theory. 

Moreover, from one such calculation, namely, a 2 x 8 ladder at a single value of A, one 
can extract the perturbative coefficients for the ground-state energy and the spin gap up to 
7th-order. That is, the ground-state energy per site is given by 

E /N =--- -A 2 - —A 3 + —A 4 + —A 5 + —A 6 - —A 7 , (24) 
0/ 8 16 32 256 512 2048 2048 ' V ; 

while for the spin gap we obtain: 

A = 1 - A + -A 2 + ^A 3 - ^A 4 - —A 5 - —A 6 + —A 7 . (25) 
2 4 8 128 1024 2048 V ; 

Of course, these two expressions can be used for arbitrary values of A, provided these 
values lie within the radius of convergence of the perturbative expansions (more on this point 
later). In Fig. || we have displayed the bulk values for the ground-state energy per site and 
the spin gap to various orders in perturbation theory; we have also included exact results 
obtained on a 2 x 8 ladder. We observe that in the strong-coupling regime (A < 0.5) the 
convergence of our results is very good. This should not come as a surprise, as the rung basis 
has been customized to handle the strong-coupling region. In contrast, the convergence of 
the perturbative results in the intermediate- to isotropic-coupling region (A <; 0.8) is harder 
to assess. However, the fact that there is no systematic improvement — nor convergence — in 
the results as one goes to higher orders in perturbation theory, strongly suggests that the 
expansion might no longer be convergent. 

In an effort to establish quantitatively the radius of convergence of the perturbative 
expansion, we have compared high-order perturbative and exact results for finite ladders 
over a range of A values. For example, Fig. [| shows the ground-state energy per site through 
50th order in perturbation theory for 4-, 6- and 8-rung ladders as a function of A. The 
radius of convergence is determined by the onset of wild excursions in these quantities at 
A c ~ 0.8. These perturbative results are indistinguishable from the Lanczos energies for 
A < A c . As the value of A c appears to have very little dependence on the size of the ladder, 
we assume it is also appropriate to the bulk system. Thus, we conclude that the radius 
of convergence of the perturbative expansion in the rung basis is A c = 0.8. This result 
indicates that perturbative expansions in the rung basis are unsuitable for dealing with the 
physically-interesting case of isotropic (A = 1) ladders. 

Finally, we conclude this section by listing in Table | the perturbative coefficients for the 
one-magnon dispersion relation, as given in Eq. (|20|). 
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B. The Plaquette Basis 



Plaquette-based perturbative results for the ground-state energy per site and for the 
spin gap appear in Table ||. By carrying out Rayleigh-Schrodinger calculations for isotropic 
(Jii = Jj_) ladders as large as 2 x 12 (i.e., six plaquettes), we have been able to determine 
the ground-state energy per site through 11th order and the spin gap through fourth order. 
Through 11th order, the ground-state energy is within 0.02% of the DMRG f|| value while 
at fourth order, the perturbative estimate for the spin gap agrees with the DMRG value to 
better than 10%. The perturbative energy per site is plotted versus the order of perturbation 
in Fig. [5] where the rapid convergence to the DMRG value is apparent. A similar plot for 
the spin gap is shown in Fig. ^. Here convergence is less dramatic presumably due to the 
relatively low order of perturbation theory to which we are restricted. As became apparent 
in the previous section, there is no point in presenting a similar table or figures for rung-basis 
perturbation theory, as it does not converge for A = 1 which is the only point where direct 
rung versus plaquette comparisons are possible. 

To extract the radius of convergence in plaquette-based perturbation theory we have 
again compared high-order perturbative and exact results for finite ladders over a range 
of A values. Indeed, Fig. [7| shows the ground-state energy per site through 50th order in 
perturbation theory for 4-, 5- and 6-plaquette systems as a function of A. Notice that for the 
present case the onset of the wild excursions is now at A c ~ 1.4; once more the agreement 
with exact calculations using direct diagonalization or Lanczos methods is essentially perfect 
for A < A c . Fig. [8] shows similar plots for the lowest triplet state of ladders with 4 and 5 
plaquettes. On the basis of the behavior shown, we estimate that A c ~ 1.25 for this state. 
It is therefore straightforward to conclude that the plaquette basis offers, to date, the only 
reliable perturbative approach for the treatment of the physically interesting || case of 
isotropic (A = 1) spin ladders. 

Finally we conclude this section by presenting in Fig. |9] the perturbative one-magnon dis- 
persion relation for the plaquette basis. The curves presented in these figures were generated 
by inserting the coefficients appearing in Table |T| into Eq. ( E?0|) . 

V. SUMMARY AND CONCLUSIONS 

We have employed rung and plaquette bases to study bulk properties of Heisenberg 
spin-1/2 ladders using Rayleigh-Schrodinger perturbation theory. In both cases, the mere 
selection of an angular momentum coupled basis enables one to separate the Hamiltonian 
into a unperturbed part, which is solved exactly, and a "residual interaction" which is treated 
in perturbation theory. In this way the starting point for the calculation are unperturbed 
wavefunctions consisting of a direct product of spin singlets for the ground state, and a 
degenerate "one-magnon" band in which one spin singlet has been turned into a triplet. 
The unperturbed ground-state energies per site in the rung and plaquette bases are given 
by —3/8 and —1/2, respectively; the spin gap starts at A = 1 in both cases. 

The one-dimensional nature of the problem enabled us to employ powerful methods for 
the extraction of high order perturbative coefficients for the bulk system from calculations on 
relatively small ladders. Indeed, using very general arguments — related to the linked-cluster 
expansion and to the scaling of observables in the bulk limit — we have shown that from a 
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calculation on a ladder with N clusters (a cluster being a rung or a plaquette) we can extract 
all perturbative coefficients for the ground-state energy and for the spin gap up to order N—l. 
Moreover, with a simple modification, namely using open-boundary conditions, the ground- 
state energy per site could be computed up to order 2N — 1. Yet, because of degeneracies 
in the one-magnon band, perturbative coefficients for the one-magnon dispersion relation 
could only be extracted up to order N — 2 by combining results from Rayleigh-Schrodinger 
calculations for ladders with N and N — l plaquettes. 

For the rung basis we have carried out calculations on a relatively small 2x8 ladder; 
this sort of computations can now be carried out routinely on a personal computer. Yet, 
the power of the method allowed us to evaluate the ground-state energy and the spin gap 
to 7th-order in perturbation theory. Perhaps as important, we could reliably estimate the 
radius of convergence of the perturbative expansion in the rung basis. This is essential, 
as such expansions are routinely assumed to be meaningful in the isotropic (J\\/J± ~ 1) 
limit We have found a radius of convergence for the perturbative expansion in the 

rung basis of only A = J\\/J± — 0.8. 

Similar calculations were also carried out in the plaquette basis. Here the parameter A 
scales the inter-plaquette to the intra-plaquette interactions; the isotropic (Jy = J±) limit 
was obtained by setting A = 1. We were able to perform calculations for the singlet ground 
state on a 2 x 12 ladder. For the triplet states, we were restricted to a 2 x 10 ladder. 
Even when a ladder of such large size was employed, the spin gap could only be computed 
up to 4th-order in perturbation theory; the ground-state energy was computed with open- 
boundary conditions so llth-order results were reported. At first sight this appears as a 
drawback relative to the rung-basis calculations described earlier. However, on close exam- 
ination one observes that rung-basis perturbation expansions in the physically-interesting 
case of the isotropic limit (A = 1) are not feasible even in principle. In contrast, the isotropic 
point is well within the radius of convergence — estimated to be at A c ~ 1.25 — for plaquette- 
based expansions. Finally, by a combination of 4- and 5-plaquette results, we were able to 
extract the one-magnon band dispersion relation through 3rd order in perturbation theory. 

In conclusion, we have formulated perturbation theory for two-leg Heisenberg spin lad- 
ders in two angular momentum coupled bases, namely the "rung" and "plaquette" bases. 
By exploiting the connectedness of perturbative energy contributions, we have been able 
to straightforwardly extract quantities appropriate to the bulk system from Rayleigh- 
Schrodinger perturbation theory for finite ladders. We find that the radius of convergence 
in the rung basis is A c ~ 0.8 while, in the plaquette basis, it is A c ~ 1.25. Clearly, the 
latter basis is strongly preferred when treating the physically interesting case of isotropic 
spin ladders which corresponds to A c = 1 in both bases. As we have stressed elsewhere ^,^3J , 
the plaquette basis is especially well-suited to the study of many aspects of spin ladders. 
Its success in the context of perturbation theory as presented here encourages us to use it 
in other investigations including a mean- field treatment of ladders (see, e.g., Ref. [|]] for a 
rung-basis approach) and in a t — J model aimed at the very interesting question of the 
dynamics of doped ladders (again, see, e.g., Ref. [Hi] for a rung-basis approach). These 
studies will be the subject of future reports. 
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FIGURES 
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2 4 Plaquette Basis 



FIG. 1. Schematic representation of the rung and plaquette bases. The solid links represent 
those interactions which are diagonal in the basis; the dashed links represent the "perturbations". 
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FIG. 2. Rung-basis perturbation theory results through order n = 10 for the ground-state 
energy per site (left) and for the spin gap (right) for systems with 4, 6 and 8 rungs using periodic 
boundary conditions. For the ground-state energy, the triangles display a 4-rung calculation with 
open boundary conditions (see text for an explanation). Here A = 0.25 and the DMRG values are 
from Ref. lO 
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FIG. 3. Ground-state energy per site (left) and spin gap (right) in the bulk limit as a function 
of A in the rung basis. Various orders of perturbation theory are displayed as well as the exact 
results on an 8-rung ladder. The DMRG values Q for the bulk system at the isotropic point 
(A = 1) are displayed with a solid square. 
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FIG. 4. Ground state energy per site for 4-, 6-, and 8-rung ladders at 50th order in perturbation 
theory as a function of A. We conclude that the radius of convergence of the rung-basis perturbative 
expansion is about A c = 0.8. 
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FIG. 5. Ground state energy per site in the bulk limit versus order of perturbation theory in 
the plaquette basis. The inset shows the same results but on an expanded scale. The DMRG value 
is from Ref. [01. 
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FIG. 7. Ground state energy per site for 4-, 5-, and 6-plaquette systems at 50th order in per- 
turbation theory as a function of A which scales the inter-plaquette interaction (see the subsection 
on Rayleigh-Schrodinger Perturbation Theory). We conclude that the radius of convergence of the 
plaquette-basis perturbative expansion is about A c = 1.4. 
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FIG. 9. One-magnon dispersion relation at various orders of plaquette-basis perturbation theory. 
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TABLES 



TABLE I. Rung-based perturbation theory for the one-magnon dispersion relation; see Eq. (20). 



Order n A^ A^ A^ Aj n) 

1 

1 1 

2 3/4 -1/4 

3 3/8 -1/4 -1/4 1/8 

4 -13/64 -5/16 -5/160 1/8 -5/64 



TABLE II. 


Plaquette-based perturbation theory for the bulk energy per 


site and spin gap of 


isotropic 2-leg 


spin ladders, [see Eq. (11)] DMRG results are from Ref. 0. 




Order n 




aw 





-0.5 


1 


1 





-2/3 


2 


-0.05815972 


0.01967593 


3 


-0.01739728 


0.06493538 


4 


-0.00329927 


0.04306155 


5 


-0.00040692 




6 


0.00021569 




7 


0.00046854 




8 


0.00041483 




9 


0.00021860 




10 


0.00005215 




11 


-0.00003610 




Total 


-0.57792948 


0.46100619 


DMRG 


-0.57804314 


0.504 



TABLE III. Plaquette-based perturbation theory for one-magnon dispersion relation; see 
Eq- ©■ 



Order n A^ A^ A^ 

1 

1 -2/3 

2 145/864 -1/9 -1/27 

3 0.10022425 0.04997348 -0.04822531 -1/27 
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